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Abstract 

The behavior of a bidisperse inelastic gas vertically shaken in a compartmentalized container 
is investigated using two different approaches: the first is a mean-field dynamical model, which 
treats the number of particles in the two compartments and the associated kinetic temperatures 
in a self-consistent fashion; the second is an event-driven numerical simulation. Both approaches 
reveal a non-stationary regime, which has no counterpart in the case of monodisperse granular 
gases. Specifically, when the mass difference between the two species exceeds a certain threshold 
the populations display a bistable behavior, with particles of each species switching back and forth 
between compartments. The reason for such an unexpected behavior is attributed to the interplay 
of kinetic energy non-equipartition due to inelasticity with the energy redistribution induced by 
collisions. The mean-field model and numerical simulation are found to agree qualitatively. 

PACS numbers: 02.50.Ey, 05.20.Dd, 81.05.Rm 
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I. INTRODUCTION 



Granular fluids are currently attracting growing interest in view of their unusual prop- 
erties, some of which are not fully understood yet. Moreover, new experiments continue to 
reveal unexpected phenomena, which have no counterpart in molecular fluids. Clustering, 
shear instability, non-Maxwellian velocity distributions, long range velocity correlations, 
non-equipartition in a binary mixture are just a few of these peculiarities Re- 
cently, another fascinating phenomenon was reported. It is the so called "Maxwell sand 
daemon" experiment ^], where a system consisting of inelastic particles enclosed in a two- 
compartment container is shaken vertically. Particles can flow from one compartment to 
the other through a small orifice located at a certain height from the basal vibrating plate. 
For strong shaking, the right and left populations are statistically equal, whereas for weak 
shaking the system spontaneously breaks the left-right symmetry. The mechanism behind 
such an unusual ordering process is the clustering induced by inelasticity. 

In fact, an imbalance in populations induced by a fluctuation can be amplified, since it 
causes a larger energy dissipation on the overpopulated side, thus suppressing the outflow 
from that compartment. At the same time, inflow from the underpopulated compartment 
is enhanced because of its lower occupation, which results in higher kinetic energies per 
particle Q|. 

Quite recently, the Twente collaboration jy] investigated the behavior of a bidisperse 
granular mixture of small and large particles using a similar experimental setup. Their 
experiments demonstrated that a bidisperse compartmentalized granular mixture has a ten- 
dency to cluster competitively. Depending on the shaking strength, one can observe different 
asymptotic configurations. 

Theoretical treatments ofgranulargases in compartmentalized systems range fro m p he- 



xgases m compartmentalized systems range irompne- 
Il2| . to molecular dynamics jl3L Il4| . to 



nomenological flux models [a, |y, 17 

more refined kinetic approaches [la]. Whereas the full solution of the inelastic Boltzmann 
equation remains a formidable task, a simple set of mean-field dynamical equations can be 



derived 



16|,[l2|. 

According to this method, the kinetic temperatures and the occupation 



numbers in each compartment are assumed to be the only relevant dynamical variables 
and treated on equal footing, a technique which naturally lends itself to capture the more 
complex phenomenology expected in mixtures. 
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Here we study a binary mixture of inelastic hard disks in a two-compartment system. 
The two species have different masses, are subjected to gravity and driven by a vibrating 
base. We extend the mean-field treatment of ref. 16] and compare its predictions with the 
results of even-driven simulations. 

The paper is organized as follows: after introducing our model in section II, we develop our 
mean-field treatment in section III. In section IV we report the predictions of our model and 
summarize our findings with a mean-field "phase diagram" , where the boundaries between 
different regimes are studied as function of the control parameters. In section V we turn 
to the event-driven simulation and find qualitative agreement with the previous picture. 
Finally, in section VI, we present our conclusions. 

II. THE MODEL 

Let us consider a two-dimensional rectangular container with horizontal and vertical sides 
of length L x and L z , respectively, divided in two equal compartments by a wall of height 
h < L z . The box contains a mixture of M\ and M2 inelastic hard disks of diameter a and 
masses mi and 1712, respectively. The disks are subjected to a gravitational force acting 
along the negative z direction and are fluidized by the sawtooth-like movement of the base, 
which oscillates with frequency v and amplitude A. The side and top walls are fixed, and 
particles collide with them in a perfectly elastic fashion. The collisions between particles are 
inelastic, and will be described by means of a velocity-dependent coefficient of restitution 



where V n is the pre-collisional relative velocity along the direction joining the centers of the 



values of relative normal velocities, the function a(V n ) assumes a constant value r, while for 
V n — > the elastic behavior a(0) = 1 is approached according to a power law. We checked 
the robustness of our results with respect to variations in t> - F° r the sake of simplicity, 
the law a(V n ) is the same for all types of collision. Binary collision between a particle of 
species fi and another particle of species v change pre-collisional velocities vi and V2 into 





two particles, v 



^fgd is a cutoff velocity and g is the gravitational acceleration. For large 
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post-collisional velocities v' x and v 2 according to: 

v[ = vi - fi Vfi (l + a)(n • vi 2 )n (1) 
V 2 = v 2 + /x /1J/ (l + a)(n- v 12 )n (2) 

where // M „ = m^jivn^ + m„), n is a unit vector directed from the center of the particle of 
type fi to the center of particle u, and vi 2 = vi — v 2 . Inter-particle collisions conserve 
momentum, but result in an energy loss proportional to (1 — a 2 ). 

The space of the control parameters is large since the system properties are functions of 
several dimensionless quantities, such as the coefficient of restitution r, the total number of 
particles M = M\ + A/" 2 , the mass ratio, the ratio of lengths L x /a, the ratio between the 
typical energy transferred from the piston to the particles and their gravitational energy, and 
so on. We study the system behavior by considering a combination of the above parameters, 



X 



whose relevance has been pointed out in |8j. The dimensionless parameter R decreases as 
the driving intensity Av increases, while increases as the dissipation becomes larger. 

III. MEAN-FIELD THEORY 

In this section we shall extend the mean-field treatment of the compartmentalized inelastic 
gas, which was introduced in refs. If], [itJ , to include the case of different species, non-equal 
occupation numbers and different partial temperatures of the system. In a realistic scenario, 
the presence of strong inhomogeneity and of the ensuing gradients make the analysis of 
the Boltzmann transport equation far too complex. The Twente collaboration proposed a 
simple flux model U while in ref. a simple coarse-grained version of the Boltzmann 
equation was introduced, leading to equations for the occupation numbers that appear to 
be very similar to those of the flux model. Such an approach has also the advantage of 
treating occupation numbers and granular temperatures on equal footing. In the following, 
we extend this method to the case of a granular mixture. 

In order to obtain our mean-field equations we shall make a series of simplifying assump- 
tions: 
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1. the system properties are regarded as homogeneous within each compartment, so that 
only the numbers of particles and the kinetic energy of the two species in each com- 
partment are the only relevant variables; 

2. the coefficient of restitution is velocity- independent, i.e. a(V n ) = r; 

3. the driving mechanism is spatially homogeneous and representable by means of a 
stochastic thermostat; 

4. the velocity distributions are Gaussian, and their variances are derived in a self- 
consistent fashion; 

5. the effect of collisions is adequately described by the inelastic homogeneous Boltzmann 
equation; 

6. particles whose kinetic energy exceeds an assigned threshold switch to the adjacent 
compartment at a fixed rate r" 1 . 

The assumption of spatial homogeneity allows us to write the phase-space distribution func- 
tion for species \i (with /i — 1, 2) as 

/ M (r,v,t) = / 0)M (v,t) . 

Here a = 1 when point r falls within the left side of the container (relative to the barrier) 
and a = 2 when r falls on the right side. Such a coarse-graining procedure allows us to write 
the following Boltzmann-like equation: 

2 

9tfaA y > f ) = XI VM/a,M> M + Bfw + *M /<»,**» fbj , (4) 
v=\ 

where is the collision term which describes the effect of inelastic collisions among particles 
belonging to the same compartment, Bf a ,n represents the action of the stochastic driving 
force associated to the heat bath, and X\y\f a ^, fb tfJ ] represents the flow of particles of species 
\x between the two compartments. 

The integral describing collisions between a particle of type \i and a particle of type 
v can be represented as 

I»Avl\fa,»Ja,A =V J dw 2 j (ffi(n • Vi 2 ) (J^faAv'OfaA^) ~ /o >A i(Vi)/a )V (v 2 )^ , (5) 



where doubly primed symbols stand for pre-collisional velocities and the r 2 factor in the 
denominator stems from inelasticity [3]. 

The number of particles of species n in compartment a is defined as 

n a ,fi{t) = J dr J dvf a>II (v,t) . (6) 

Similarly, we define four partial granular temperatures: 

/f Tfl V 2 
dr / dw— ^— f a ,p(v,t) . (7) 

Now we derive the governing equations for occupation numbers and kinetic temperatures. 
By integrating eq. (jlj) with respect to velocity v 1; and recalling that both the collision 
integral and the heat bath term conserve the number of particles in each compartment, one 
can see that the rate of change of n a At is determined solely by the exchange term X: 

d t n a ^(t) = J dvX[v\f a ^,f b J ■ (8) 

According to the assumption of point 6 (see list above), the exchange term has the form 

XWa^ftJ = --9(\V\-U s )[f a ^,t) - / M (v,t)] . (9) 
Ts 

In other words, one assumes that particles having velocities larger than a given threshold u s 
cross the barrier with probability t^ 1 per unit time. The resulting equation for n a ^ is 

^^ = -^/ dv [/^( v ^)-/^( v ^)Ml v l-^)- ( 10 ) 

Mathematical convenience suggests an additional assumption for the form of the heat bath 
operator B: we replace the effect of the periodic vibration of the piston with a sequence of 
uncorrelated random kicks whose amplitude is distributed according to a Gaussian j^, I2]] . 
The effect of this stochastic acceleration is described by the Langevin equation 

dvj , , 

where £ is a a Gaussian white noise with zero average, and variance given by 

Ut)Ut')) = v%At-t'). (12) 

Here if measures the intensity of the stochastic driving. This leads to the following forcing 
term: 

rf d 2 

Bf a ^ = YQ^faAV't) ■ ( 13 ) 
6 



Notice that particles receive an energy proportional to their masses. Next, we consider the 
time evolution of partial kinetic energies T^, which is obtained by multiplying eq.(jU by 
m fl v 2 /2 and integrating over velocity. 



yd t {n a]tl T aitl ) = ~y22 / c?vv%„[v|/ aiM , f aiU ] + -± j dwv 2 Bf a ^ 

-^/dvw 2 [/a rf .(v,t)-/Uv,*MM-tt.), (14) 

Finally, in order to derive an explicit expression for the governing equations we assume 
Gaussian velocity distributions: 

(15) 

Inserting eq. (|T5|) into equation (fTUj) we obtain 

= - [n b ,,e- T ^ - n a ,,e- T ^] , (16) 
at t s 

where the temperature T^ s is given by T^ s = ^m^ul and the right-hand side represents the 
difference between the incoming flux and the outgoing flux of compartment a, as in the flux 
model. 

The calculations involving the collision term are quite lengthy, but straightforward. Here 

231 ] . The final form 
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we will not report the detailed derivation, which can be found in refs. 
of the governing equation for the granular temperature (say, of species 1 in compartment a) 
is: 

dt r s 

- - [(n a ,ie- 2 W T - 1 - n„,ie- Tls/T M)(2T ls - T a>1 

rr ( T a ,l T a 2 

- 7ll«o,l-to,l - 7i2^in ,i 1 

\ mi ^2 

T — T 

- Ki 2 min a ,i( ■ ) +771x7] n a ,i , (17) 

mi + 7TL2 

The corresponding equations for other choices of species and compartment are obtained by 
simple substitution of indices, where the coefficients •y fJiU and are given by 
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^ = 2(l + r)a^— ^ + . (19) 

The constant A is the area available to the particles in each compartment. 7^ takes into 
account the effect of collisions between particles belonging to species fi and u, respectively. 
Interestingly, does not vanish in the elastic limit (r — > 1), since it is associated with 
the restoring mechanism which tends to equalize the partial kinetic temperatures of two 
different components in a standard fluid. 

The complete set of equations (|TT)|) and (jl7j) consists of eight non-linear coupled equa- 
tions, of which only six are independent, because the occupation numbers - being globally 
conserved - have to satisfy ni tfl + = M hl . 



IV. THE MEAN-FIELD SCENARIO 



To be able to compare the predictions of the proposed mean-field model against the results 
of numerical simulation, we need to introduce a suitable definition of R (eq. EJ) for the case 
of our thermal bath. Specifically, we rewrite the first factor of equation © as mgh/(P 2 /m), 
where P is the typical momentum transferred to a particle that collides with the piston. 
We replace the energy scale mgh - associated to barrier crossing - with the analogous scale 
T^s, and the transferred momentum scale P with its equivalent for the stochastic forcing of 
eq. El Tnrj^/f^. Thus, for the mean-field model, the analogous of eq. El can be written as: 

T N 2 a 2 
RMF = Z ^{l-r 2 Y-^ r . (20) 



X 



Obviously we expect the correspondence between R and Rmf to be useful only for qualitative 
comparison of the mean-field model against numerical simulation. 

We shall now consider some of the properties of the asymptotic solutions of eqs. (jl6|) and 
(I17|) . The specific numerical values of the control parameters correspond approximately to 
those used in the experiments by (|. 

A symmetric solution, characterized by equal state variables in the two compartments, 
n\,n = and Ti jAt = T 2jAt , always represents a fixed point for the governing equations. 
However, as the parameter Rmf increases, such a symmetric fixed point becomes unstable. 
At small driving intensity or high dissipation the solution ceases to be symmetric and the 
particles of both species tend to cluster in one of the two compartments. Such a symmetry 
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breaking behavior was studied experimentally and explained in terms of flux models 
By numerically integrating eqs. ()16|) and (fTTj) . we are also able to observe the same simmetry- 
breaking behavior for the case of our mean-field description. 

In fig. ^ we show the temporal evolution of occupation numbers for the two species (in a 
single compartment) for different choices of the driving intensity. We observe two dynamical 
scenarios. In the first scenario, the occupation numbers approach a stationary symmetry- 
broken solution via a simple bifurcation (i.e a critical point). This corresponds to the 
situation considered in ref. 2^ |. and occurs for small mass differences or small concentrations 
(m^jnix = 2, for example). In the second scenario, the temporal evolution of the occupation 
numbers approaches a limit cycle so that that particles of each species oscillate back and 
forth between compartments. This occurs when the difference in mass is large, or when 
masses are equal but sizes are different enough, in analogy to the observations of Lambiotte 



et al. 



24 ]. 



Although the system is far from equilibrium and has a finite extension we shall employ 
the term "phase diagram" in order to stress the existence of different dynamic regimes. The 
resulting phase diagram in the plane Rmf vs. 1112/ mi is shown in fig. El No oscillations can 
be observed for small mass asymmetry ~ 1) and the crossover from the symmetric phase 
to the asymmetric phase is similar to what occurs in the case of a one-component system. 
However, when the mass asymmetry increases (m 2 /mi > 3) an intermediate oscillatory 
regime appears, and no stationary solution is attained anymore. 

The effect of concentration on the appearance of oscillations is instead shown in fig. 01 
For small concentrations of the heavy species A2/ [M\ + A/2) there is a direct crossover, as R 
decreases, from the asymmetric "phase" to the symmetric "phase" . When the concentration 
A2/(A/i + A/2) of the heavy particles increases, there appears an island of the oscillating 
"phase" , which subsequently disappears when the concentration of heavy particles becomes 
too large. The oscillation period decreases as Rmf decreases and is larger near the boundary 
between the "broken phase" and the "oscillatory phase" . In the following, in order to simplify 
the notation, we shall denote N\ = n\ \ and N% = ni\- 
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FIG. 1: Four different behaviors obtained by varying the dimensionless parameter Rmf, and 
keeping the mass ratio m^jvnx = 8 and the population ratio N1/N2 = 70/20. The continuous 
line (iVi) refers to the light species and the dotted line (N2) to the heavy species, both in the 
right compartment. In panel A (Rmf = 6.0) the asymptotic solution is symmetric. In panel B 
(Rmf = 12) the occupation numbers oscillate in time. Panel C (Rmf = 24) shows a case where the 
solution displays oscillations with longer periods. Finally, panel D (Rmf = 48) illustrates a typical 
simmetry breaking solution, where occupation numbers in the two compartments are different, for 
each species. The time is measured in units t s . 

V. NUMERICAL EXPERIMENTS 



In order to verify whether the hard disk model correctly describes experimental observa- 
tions [?J and also to probe the qualitative validity of our mean-field scenario, we now turn 
to discuss a set of numerical experiments. We confirm that the behavior of the mixture, 
as already revealed by mean-field calculations, strongly depends on the specific choice of 
control parameters. 

The simulated dynamics consists of a succession of ballistic trajectories and inter-particle 
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m 2 /m i 



FIG. 2: Phase diagram for a system with M\ = 70 light particles and A/2 = 20 heavy particles. 
The transition values of Rmf are plotted as functions of the mass ratio m,2i 'm\. 

or wall-particle collisions. With respect to the code employed previously (see references jl^ 
and [26]) we introduce a dividing, elastic barrier, visibile in fig. |3J The two species are chosen 
to be smooth rigid disks of equal diameters a = 0.5 cm and unequal masses mi = 6.545-10~ 2 g 
and m2 = 5.236- lO -1 ^, respectively. The mixture is composed of A/"i = 60 light particles and 
A/2 = 30 heavy particles. The vibrating base is driven according to a symmetric sawtooth 
waveform and the vibration amplitude is set to A = OAa. Finally, the width L x of the 
compartment ranges from 34cr to 240cr, while its vertical dimension is L z = 82.4(7, and the 
height of the dividing wall is h — 12a. 

Our observables are computed as instantaneous averages in the case of global quantities, 
such as the occupation numbers for the two species and the partial kinetic energy per particle, 
and as time averages (ergodic averages) for other quantities, such as densities, temperatures 
and velocity distribution functions. 

For strong driving (i.e. small R, as already shown in the mean-field case) the system 
approaches a stationary state where the number of particles in each compartment, or average, 
is the same. Such a steady state is maintained by a continuous exchange of particles of both 
species between the compartments. The average kinetic energies are sufficiently high so 
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FIG. 3: Phase diagram for a system with Mi + Mi = 90 particles as function of M2/ (Mi + A/2), for 
a mass ratio mi/ mi = 4. 
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FIG. 4: Typical snapshots of the process described in Section V. Open circles indicate the light 
particles, black circles the heavy particles. System parameters are v = 47 Hz, L x = 34a, mi/ mi = 
8, Mi = 60, A/2 = 30, A = 0.4<7 and <r = 0.5cm. 

that particles can jump over the barrier. Upon increasing R the symmetry is spontaneously 
broken: the system can either attain an asymmetric stationary state, or approach a limit 
cycle where populations keep switching back and forth between compartments in a coherent 
fashion. These scenarios are shown in figures El for mi/mi = 8, r = 0.85 and v = 47 Hz, 
where the crossover from one regime to the other occurs by varying the value of L x . For 
L x = 160<t and L x = 120a the particles are nearly equidistributed among the compartments. 
On the other hand, when the system is narrow enough (L x = 62a) the particles - after an 
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initial transient - tend to localize in the same compartment for time intervals longer than our 
observation time. However, when we change the width of the compartment from L x = 62<j 
to L x = 8O0", pronounced oscillations appear and the populations appear to move back and 
forth. The global picture is consistent with our previous mean-field analysis. 

We also explore the effect of varying the driving frequency at fixed L x , and we find that 
the breakdown of the symmetric stationary state still occurs, as shown in fig. |U1 Moreover, 
the oscillatory behavior only appears for large mass asymmetry. On the other hand, the 
time series of fig. [7| shows that the corresponding mass asymmetry (m 2 /mi = 2) is not 
large enough to induce clean oscillations, so that only the phenomenology of one-component 
driven granular gases can be observed. 

Finally, we remark that the same kind of oscillations may occur when the diameters of 
the two species are different enough, even in the case of equal masses (data not shown). 

Let us consider in more detail the oscillatory behavior we observe (for example, refer to 
panel C of figure EJ). After an initial clustering of both species in a compartment - say the 
left - a net rightward flux of light particles establishes and persists until a sufficient number 
of them have changed compartment. At this point the heavy particles, too, start jumping 
to the right, eventually creating in the right compartment a cluster of both species which 
is totally similar to the initial situation of the left compartment. After reaching this stage, 
the process repeats itself in the opposite direction. 

Such a peculiar behavior has its origin in the mass and/or size difference between the 
two species and in the associated breakdown of energy equi partition that occurs in a binary 



vibrated granular gas [23, |2Sl l29j |. Krouskop and Talbot |3j| studied how the mass ratio 
affects the exchange of energy of each species and computed the average energy change of 
a particle of species fi colliding with a particle of species v. Assuming that each component 
has a Maxwell velocity distribution and a partial temperature T M , the average energy change 
in 2D is 

(21) 



AP _ a + r) 



I _|_ r \ ^ n — _ 2—^- 

m^rrif, + m u ) 



2V2 

Such an equation shows that for T 2 /Ti > 1 a light particle - on average - gains energy on 
colliding with a heavy particle. In our case, the imbalance in temperatures is guaranteed 
by the fact that heavy particles receive more energy than light ones when colliding with the 
vibrating base. Such an imbalance is roughly proportional to the mass ratio. Thus heavy 
particles are able to gain more energy from the driving base, and they can also tranfer it 
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FIG. 5: Number of light particles (solid line, N\) and heavy particles (dotted line, iVjj) in the left 
compartment, as a function of time, for different values of L x : L x = 160a (A), L x = \20a (B), 
L x = 80<T (C) and L x = 62a (D). System parameters: r = 0.85, v = 47 Hz (period T = 
A = 0.4a, a = 0.5cm and m 2 /mi = 8. The panels correspond to i? = 5.8 (A), 10.4 (B), 23.4 (C), 
39.0 (D). 

to light particles, effectively heating them. Since the system is not in equilibrium, kinetic 
temperatures depend on the particle numbers in each compartment and are not known a 
priori. 

In order to check the validity of the above argument, we focused on the steady-state 
dynamics between transitions, and numerically evaluated the partial temperatures T\ and 
T2 in the compartment hosting the majority of heavy particles. Such temperatures mainly 
depend on the number iV of light particles, while they are not signifcantly affected by 
the number of heavy particles, because the latter population is basically constant between 
transitions (see fig. EJ). On inserting the above partial temperatures into equation (j2H), and 
weighting the energy changes with the numerically evaluated collision frequencies, we obtain 
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FIG. 6: Number of light particles (solid line, A^) and heavy particles (dotted line, A^) in the left 
compartment, as a function of time, for different driving frequencies v (period T = 1/u): v = 37 
Hz (A), v = 47 Hz (B), v = 50 Hz (C), v = 52 Hz (D). System parameters: r = 0.85, L x = 68a, 
A = 0.4ct, a = 0.5 cm and m\/m2 = 8. 



the average energy changes per particle and per collision, separately for each species: 

AEnMu + AE 12 M 12 



and 



AK 1 



M u + M12 
AE 21 M 12 + AE 22 M 22 



(22) 



(23) 



M 22 + M 12 

Here M M!/ represents the number of collisions per unit time between species \i and v. Fig. El 
shows AKi and Aif 2 as a function of A/". We observe that light particles, on average, gain 
energy. Such a gain increases as N decreases, because in this situation M\ 2 is much larger 
than M\\. Thus light particles are more likely to jump into the empty compartment, leaving 
the heavy ones behind. At this point the heavy particles - no longer transferring energy to 
the other species - become more energetic and also move to the other compartment, joining 
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FIG. 7: Number of light particles (solid line, N\) and heavy particles (dotted line, N2) particle in 
the left compartment, as a function of time, for different choices of L x : L x = 240a (A), L x = 220cr 
(B), L x = 200cr (C) and L x = 160a (D). System parameters: v = 47 Hz, A = 0.4<r, a = 0.5cm and 
m 2 /m 1 = 2. The panels correspond to R = 2.6 (A), 3.1 (B), 3.7 (C), 5.8 (D). 

the light particles once more. 

Let us now consider in more detail the statistical properties of the system by analyzing 
a particular case shown in fig. |HJA.. One can observe that the occupation numbers fluctuate 
around well defined plateau values for periods as long as several thousands driving cycles 
(of duration T = 1 jv) . 

Interestingly, the distribution of the horizontal position of the center of mass displays a 
double peak structure, indicating that the system spends most of the time in configurations 
where the majority of the particles is on one side - an asymmetric situation. 

We can better study the switching dynamics of figJBl^- by defining a collective residence 
time r as the period during which at least 80% of the heavy particles persist in a given 
compartment. The ensuing distribution of residence times is shown in fig. The average 
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FIG. 8: Panel A: time evolution of the occupation numbers N\ and N2 in the left compartment. 
Panel B: probability distribution of the horizontal position x of the center of mass (whole system). 
The left-right asymmetry is due to the limited size of the sample used to compute the histogram. 
System parameters are v = 47 Hz, A = 0.8a, L x = 68a and tt^/toi = 8. 



n order 



2J and 



residence time is about 6 • 10 3 T and the distribution appears to be exponential, 
to ascertain whether the oscillations are periodic, as reported by Lambiotte et al. 
predicted by our mean-field model, we study the signal of fig. |H] in the frequency domain. 
The resulting power spectrum P(uj) is shown in the inset of fig. El and displays no evidence 
of periodicity. As a further check, we generated a "synthetic" bistable signal by assuming 
that the occupation numbers of the heavy species in each compartment can assume only 
two values (say N 2 and 0), and that the residence times are randomly distributed according 
to an exponential density distribution whose characteristic time is the same as the one 
we measured. The corresponding power spectrum is also shown in inset of fig. and 
appears to be very close to the original one. Thus the numerically observed bistability is not 
periodic, but has a stochastic nature. Such a stochastic switching resembles what occurs in 
thermally activated processes Our mean- field model is of course unable to reproduce 
non-deterministic features, but it can still account for the non-stationary, bistable character 
of the reported dynamics. 

Since the fast particles play a major role in activated processes, we computed the veloc- 
ity probability distributions for the two species, in each compartment. Such distributions 
are shown in fig. El for the horizontal component of the velocity. On rescaling by the 



17 



I 




20 40 60 

N 

FIG. 9: Average energy change for particle-particle collisions, as a function of the number of light 
particles ./V present in the compartment where the majority of the heavy particles is located. AKi 
(open circles) and AK2 (squares) correspond to the energy change for light and heavy particles 
respectively, in the full compartment. 

corresponding mean squared velocities, they appear to collapse on two curves, one for the 
"empty" compartment and another for the "full" one. The two distributions deviate from a 
Gaussian law and can be fitted by f(v x ) = A/ (exp(/3v^) + exp(—/3v^)) (as already found in 
ref with the following values for the fitting parameters: Af = 0.9, /3/ = 1.38, A/ = 1.14 
for both components in the "full" compartment and Ae = 0.8, fie = 1.1, Ae — 1-35 and 

(2) (2) (2) 

A e = 1.03, pi. = 1.6, Ae = 0.98 for the component 1 and 2, respectively, in the "empty" 
one. It is interesting to observe that the exponent A increases with the occupation number, 
in agreement with the idea that high velocity tails are associated with those particle that 
underwent a smaller number of inter-particle collisions. 

We now turn to consider the vertical temperature and density profiles (see fig. EJ)- The 
temperature profiles indicate that the kinetic temperature of the heavy species is larger than 
that of the light one (being the values of these two quantities dominated by the bulk of the 
distributions). Remarkably, in the upper region the light particles are more energetic than 
the heavy ones. 

The density profiles show peaks at different positions for the two species, indicating that 
the light particles tend to stay above the heavy ones. One can also observe that the upper 



18 




T/T 

FIG. 10: Probability distribution of the residence time r for the (hardened) occupation number 
N2(t) of heavy particles. The dashed line is an exponential fit. In inset: power spectrum P(uj) 
of the hardened occupation number (solid line) and power spectrum of the synthetic dichotomic 
signal used for comparison in the text (dashed line). Simulation parameters are the same as in 

fig- El 

region, well above the top of the dividing barrier, contains mostly light grains. 

For the sake of comparison, we also performed simulations of a reference system with a 
slngl e co mp _ S For s-naH v alu e s „ t , tbe oW y p rofll e S ni W a nd „,(,) reSult v ery 
similar to the case of the two-compartment system. For higher values of z, both n\{z) and 
nziz) are larger than the corresponding profiles in the reference system. The local granular 
temperatures Ti(z) and T 2 (z) display a similar behavior. For small values of z, T 2 (z) is 
larger that Ti(z), because heavy grains gain a larger amount of energy when colliding with 
the vibrating base, but for large values of z this relation is reversed: light particles are 
"hotter" and their temperature profile is characterized by a slower exponential decay. This 
is consistent with the mechanisms discussed in section V. 
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FIG. 11: The rescaled velocity distribution functions v x ms P{v x ) = f(v x /v. 
ties. The open and closed symbols correspond to the "empty" and "full" compartment, respectively. 
Circles and the squares correspond to the light and heavy species, respectively. The independent 
variable v x was rescaled by its mean squared value v r x ms : 38.29cm/s (component 1) and 22.12cm/s 
(component 2) for the "full" compartment; 59.44cm/s (component 1) and 32. 40cm/ s (component 
2) for the "empty" compartment. The solid and dash-dotted lines are the fitting laws discussed in 
the text, while the dotted curve is a Gaussian distribution, as an aid for the eye. 



VI. CONCLUSIONS 



In this paper we studied the behavior of a two-dimensional, driven, inelastic granular 
mixture in a compartmentalized container. Our work is hinged on two complementary ap- 
proaches. The first approach consists in reducing the dynamics, by means of an appropriate 
coarse-graining procedure, to a simple set of non-linear coupled ordinary differential equa- 
tions for the relevant observables, namely the number of particles and the average kinetic 
energies in each compartment. These mean-field equations are solved numerically and show 
the existence of three qualitatively different regimes: 1) for "low" values of the dimension- 
less parameter Rmf the asymptotic state of the system is symmetric, i.e. the particles 
of each species are equidistributed in the two compartments; 2) for intermediate values of 
Rmf the system attains a limit cycle where the populations in each compartment exhibit 
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FIG. 12: The granular temperature (a) and density (b) profiles versus height, for the "full" com- 
partment. Square and circles correspond to light and heavy components, respectively, for the 
case with two compartments. The dashed and the dotted lines correspond to light and heavy 
components in the case of a single compartment. 

an oscillatory behavior, and segregation of species is observed; 3) for high values of Rmf 
the system approaches an asymptotically steady state, with unequal populations in the two 
compartments and no segregation of species. 

The second approach is a more realistic event-driven simulation with physical parameters 
similar to those employed in laboratory experiments. Our simulations show qualitatively 
similar scenarios. However, whereas in the mean-field description we find two stationary 
regimes, a symmetric configuration and an asymmetric one, separated by an oscillatory 
regime, our numerics indicate that the intermediate oscillatory "phase" is indeed a stochastic, 
bistable regime. In both cases, the bistability emerges only when the masses or diameters 
of the species are different enough. In addition to that, the typical time scale characterizing 
transitions in the bistable regime grows as one departs from the asymmetric phase and moves 
towards the symmetric phase. The mechanism that sustains such a non- stationary dynamics 
was explained in terms of the microscopic kinetics of a mixture of inelastic particles. 
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